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Optical Density Using a Target 


without dust 


with dust 





Average particle Size from Optical Density 


8m 
(D} = 2v ,p, In((U, -1)) ML, -ly)) 


(D) = average scattering cross-section diameter 
L =distance laser travels through dust ~ 0.5 in = 0.0127 m 
m = mass transfer rate ~ 0.0125 kg/s 
/, =luminosity of target with laser without dust 
[/ =luminosity of target with laser with dust 
[, =luminosity of target without laser (background light only) 
P, = particle bulk density = 3100 kg/m? 
v, =eV, = particle velocity m/s (v, = gas velocity = 17.7 m/s) 


| =0.299R+ 0.587 G+0.114B , luminosity according to PhotoShop 


Luminosity Measurements 


white target 


gray target 





Average particle Size from Optical Density 
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(particle velocity)/(gas velocity) 


Optical Density 2D Calibration 


Test Transparency 


wuiiuascemecemmccemeceeescees — _D = 0.050 [in] 


suai D= 0.025 [in] 


\ xCalc v2.nb 


c= ((pLp- pO) / (pL - pO)) 
xX[plp , pt_, pl _] = -Log[c] 
x[42.23, 2, 59.76] 
Exp[-x[41.9, 2, 59.02] ] 
x[(22.5, 2, 42.61] 
Exp[-x[22.5, 2, 42.61] ] 
x[11.7, 2, 19.4] 
Exp[-x[11.7, 2, 19.4]] 

-pO + pLp 

-po+pL 


-pO + | 
-pO+pL 


-Log 


0.361684 
0.699754 
0.683589 
0.504802 
0.584344 


0.557471 


Quit] 
1-Pi/és//4 
0.607301 
a=10 

L=i1 
M=10000000 


x=L/M 


Exp[-aL] // N 
(l-ax)*M // 


10 
1 


1o000000 


1 
10000000 


0,0000453999 


0.0000453997 


Limit[(1-aL/m)*m, m+ Infinity] 


en 





2D Test Image 





Optical Density Using Side-Scatter 







Luminosity Values as a Function of x 
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Luminosity, L 


90 
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to 


70 


Luminosity Values as a Function of x 


—B| x+,/(D,/2)* -(D, /2—xsin0d)” 
L(x) = Ae 


© Measured Luminosity from PhotoShop 


— Model Fit 
2 A =32.3 
B =33.8 m! 
D, =0.0191 m 
Q =20° 
Ly) = 67.393 
* 
(2 
O 
20 40 60 80 100x1 Oo” 


xX [m] 


Average Particle Size from Optical Density 
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D =average scattering cross-section diameter of particles 
B  =from luminosity fit = 33.8 m’! 

M =mass transfer rate ~ 0.0250 kg/s 

D, = diameter of spray ~ 0.75 in = 0.0191 m 

P, = particle bulk density = 3100 kg/m? 

v,=¢¥V, = particle velocity m/s (v, = gas velocity = 31.1 m/s) 


Average Particle Size from Optical Density 
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40.3 um 


m. 


(particle velocity)/(gas velocity) 


Optical Density Using Side-Scatter with 
Camera at Arbitrary Angle 





Luminosity Values as a Function of x 
12-16-11 experiment 
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R= Scwssscosy) Gsiae42) e 05 





Ronee 





oy 
oeCe ] = {€Cos[o] Cos(e], ¢Sin[d] Cos[e] -2/Tan[e], Sin[y] +=}; = 
MatrixForm[ p= [<¢]] 
HEP = {{1, 0, OF}, £0, 0, OF, (0, O, L}}; 
MatrixForm[xzF ] 
SOlve[R“2 -- (x2P.pe[(¢]).xzP.pe(¢]}, ¢] 
follz ,@, #1, 92: , A, #I}= Es. %: 
FullSimplify[pi[z, o, #], {R > 0 == 0, Im[z] =: 0, Im[@] =: 0, Im[#] =: 04] 
FullSimplify[p?[z, o, #], {R= 0 =: 0, Im[z] =: 0, Im[@] =: 0, Im[#] =: 04] 
FullSimplify[p2[z, 0, 0], {RB + 0, Im[z] =: 0, Im[@] =: 0, Im[#] =: 04] 
FullSimplify[p2[0, 0, #], {R= 0 -- 0, Im[z] -- 0, Im[@] -- 0, Im[#] =: 04] 
FullSimplify[p2[z, 0, #], {R= 0 == 0, Im[z] =: 0, Im[@] == 0, Im[#] =: 04] 
FullSimplify[p2[z, 0, Pi/2], {R= 0 =: 0, Im[z] =: 0, Im[@] -- 0, Im[#] =: 03] 
FullSimplify[p2[z, o, 0], {R= 0 == 0, Im[z] =: 0, Im[@] =: 0, Im[#] =: 04] 
‘ € Cos[p@] Cos[¥] 
| -£ Cot[é] +£Cos[#] Sin[¢] 
m4 6 Sin[4] 
‘loo; 
c o 0 | 
ool, 
{les [-16 2 Sin[¥] -.f{-4(sR'-s2") (-6-2Cos[2@] -Cos[2¢-2 4] +2Cos[2 4] -Cos[2¢+24]) +2562" sin[y]"}] / - 
(2 (64+2Cos[2¢] +Cos[2¢-2 4] -2Cos[2 8] +Cos[2e+2g]))}, 
fen [-16 2 Sin[¥] + {-4 (sR? -82") (-6-2Cos[2@] -Cos[2¢-2 4] +2Cos[2 4] -Cos[26+24]) +2562" Sinty]"}] / 
(2 (64+2Cos[2¢] +Cos[2¢-2 4] -2Cos[2 $] +Cos[2p+29])) bt 

2 [2 2 sinty] +af (R-2) (Re) (34+ Cos[26] -2Cos[2 4] Sin[g]*) +42* Sint]? | | 
- 34Cos[2¢] -2Cos[2 4] Sin[d]* 
Z [-22 sinc +af (R-2) (Rez) [34+ Cos[26] -2Cos[2 4] Sin[@]*) +42* Sin[y]* | 

DEX, ~, y) = ; it 
34Cos[2¢] -2Cos[2$] Sin[@] 

R j 
-gSin[v] + Ri- 242! sintg]! 
R-=& | 
pi -z? ahs[Cos[#]] Sec(#]! | 
Quit [] ] 


100% = 





2\-2zsinyw + ((R —z)(R+ 2)(3 +cos 2¢—2cos 2y sin’ b)+ 42° sin” vy) 


Z,9,0) = 
PEP) 3+cos 2¢—2cos 2y sin’ 


z=R-xsin@ 


L(x) = Ae Bate ew) 





Luminosity, L 


Luminosity Values as a Function of x 


© Measured Luminosity from PhotoShop 
— Model Fit 


L(x) = Ae BOte ow) 4 Ly 


= 16.8 

= 38.9 m'! 
= (0.0127 m 
= 45° 

= 120° 

Y =60° 

L, =12 
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Average Particle Size from Optical Density 
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= average scattering cross-section diameter of particles 
= from luminosity fit = 38.9 m:! 


= mass transfer rate ~ 0.0370 kg/s 
= diameter of spray ~ 1.0 in = 0.0254 m 


= particle bulk density ~ 3100 kg/m” 


v,=¢¥V, = particle velocity m/s (v, = gas velocity 31.1 m/s)" 


“ gas flow = 4.7 CFM through a 3/8 inch pipe 


Average Particle Size from Optical Density 
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29.1 um 


(particle velocity)/(gas velocity) 


Regolith Test Bin (RTB) 
Swampworks, 
Kennedy Space Center, FL 


Planetary Regolith Test Bed 





3x3 Laser 


White Target 


Reflectors 
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Intensity of Reflected Beam: I(x) = Ibe 7 @eax/4 





Extinction Coefficient: a= : In (=) 
TQox lo 


Q. = Extinction Efficiency ( ~ 2) 


x = Total Optical Distance 
= 2L (where L is bin width ~ 7.9 m) 





a = —-—— In (=) [Mm] 


lo 


«=| D*N(D)dD 
0 


D = Particle Diameter [m] 


N(D) = Particle Size Distribution [m~” m™] 


alpha(2L) [1/m] 





0.14 4. ce 
a= In | —— 
TQ 0X 
0.12 
646 = X= 2 
L = 26 ft 
0.08 Q, = 
0.06 
0.04 
0.02 


09-25-15 Test 





time [s] 


y [m] 
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x [m] 


a [m+] 
0.2 


0.15 


y [m] 


t= 12s] 





x [m] 


a [m+] 
0.2 


0.15 


y [m] 


t= 14 [s] 





x [m|] 


|| | 
+++ 


a [m+] 
0.2 


0.15 


y [m] 


t= 16 [s] 





x [m] 


et td 
|| {| | | | | | r 


a [m+] 
0.2 


0.15 


y [m] 


t= 18 [s] 
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t= 22 [s] 
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0.2 


0.15 


y [m] 


t = 26 [s] 


| 
al 
| 
aha 
= 
| 
ii 
| 


= =P 
fii} ||} ii] | lg 





x [m] 


a [m=] 
0.2 


0.15 


y [m] 


t = 28 [s| 





x [m] 


a [m+] 


0.2 


0.15 


y [m] 





t = 30 [s] 
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t= 45 [s] 
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Laser largets 


Photometers 





Left Camera — Right Camera 


Variable angle eight laser pointer plus camera (VA8LP+C). 
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Diffuse Reflection 


Retro Reflection 


Specular Reflection 


0.=Oton 





Diffuse (Gray) Target Diffuse (White) Target 


Diffuse (White) 
Diffuse (Black) Target Retro Target 





Camera image, DSC_0037.jpg corresponding to above Figures c and d: 
ISO = 100, shutter speed = 1/320 s, and aperture = f/11. 





DSC_0037.JPG 
No Filter 


Max [= 255. 
Avg [= 218.847 
Avg D = 23.4529 
Sat Cnt = 87 


DSC_0038.JPG 
No Filter 
Max I= 255. 


Avg [= 220.741 
Avg D = 23.1524 
Sat Cnt = 76 


DSC_0039.JPG 
YG Filter 


Max I= 255. 
Avg [= 217.703 
Avg D = 20.7143 
Sat Cnt = 63 





DSC_0040.JPG 
YG Filter 


Max I = 255. 
Avg [= 216.244 
Avg D = 20.0267 
Sat Cnt = 42 


DSC_0041.JPG 
2xYG Filter 
Max I = 255. 


Avg I= 197.824 
Avg D = 16.7746 
Sat Cnt = 16 


DSC_0042.JPG 
2xYG Filter 


Max I = 255. 
Avg [= 198.391 
Avg D = 16.5453 
Sat Cnt = 14 





DSC_0043.JPG 
SGI Filter 


Max I = 255. 
Avg I = 200.598 
Avg D = 16.5068 
Sat Cnt = 12 


DSC_0044.JPG 
SGI Filter 
Max I = 255. 


Avg I= 198.065 
Avg D = 16.5837 
Sat Cnt = 12 


DSC_0045.JPG 
SG2 Filter 


Max I = 255. 
Avg I= 180.729 
Avg D = 13.3512 
Sat Cnt = 1 





DSC_0046.JPG 
SG2 Filter 


Max [= 254. 
Avg I= 176.553 
Avg D = 14.3175 
Sat Cnt = 0 


DSC_0047.JPG 
SG2+YG Filter 
Max I = 197. 


Avg I= 136.62 
Avg D = 9.50789 
Sat Cnt = 0 


DSC_0048.JPG 
SG2+YG Filter 


Max [= 188. 
Avg I= 128.167 
Avg D = 10.3418 
Sat Cnt = 0 





DSC_0049.JPG 
SG2+2xYG Filter 


Max [= 130. 
Avg I= 86.6522 
Avg D = 7.65304 
Sat Cnt = 0 


DSC_0050.JPG 
SG2+2xYG Filter 
Max I= 119. 


Avg I= 82.1475 
Avg D = 8.81292 
Sat Cnt = 0 


DSC_0051.JPG 


SG1+SG2+2xYG 
Filter 


Max [= 16. 

Avg I= 11.8864 
Avg D = 7.48482 
Sat Cnt = 0 





DSC_0052.JPG 


SG14+SG2+2xYG 
Filter 


Max [ = 22. 

Avg I= 14.8148 
Avg D = 5.86323 
Sat Cnt = 0 





Image processing results for f/11 test. 























Max Saturated Spot Image | Image spot y 
Image File Tatler Pixel Count Diamete | Width | Height (Kolo) ge 
hay Fraction r [pixels] | [pixels | [pixels 
] ] 

255. 87. 0.201389 218.847 23.4529 3872 2592 2465 2058 
255. 76 0.180523 220.741 23.1524 3872 2592 2465 2059 
255. 63 0.186944 217.703 20.7143 3872 2592 2456 2056 
255. 42 0.133333 216.244 20.0267 3872 2592 2457 #2058 
255. 16 0.0723982 197.824 16.7746 3872 2592 2470 ~~ 2072 
255. 14 0.0651163 198.391 16.5453 3872 2592 2470 #2074 
255. 12 0.0560748 200.598 16.5068 3872 2592 2451 2050 
255. 12 0.0555556 198.065 16.5837 3872 2592 2452 2049 
255. 1 0.00714286 180.729 13.3512 3872 2592 2443 2047 
254. 0 0. 176.553 14.3175 3872 2592 2445 2049 
197. 0 0. 136.62 9.50789 3872 2592 2491 2063 
188. 0 0. 128.167 10.3418 3872 2592 2493 2072 
130. 0 0. 86.6522 7.65304 3872 2592 2494 #2074 
119. 0 0. 82.1475 8.81292 3872 2592 2494 2081 
16. 0 0. 11.8864 7.48482 3872 2592 2484 2098 
yee 0 0. 14.8148 5.86323 3872 2592 2485 2098 


Normalized Power, I(p) 


© data at f/36 
data at f/11 
@ datax 10 at f/11 
—— fit to f/36 data + data x10 at f/11 
- - - above fit x 0.1 


0.1 SG2 © 


ISO = 100 
ss = 1/320 
f/stop = f/36 
I(p) =A*(p-pO)“b 
pO = -127.1 ® 
A =3.5/79e-11 
b =4.048 





0.01 s 


_-7— fistop = f/11 
SG1+SG2+2YG @ _ -~ A =3.579e-12 


—_ 


_— 
—_ 
-_ — 
—_—=— — 
—_= = 


0.001 


Pixel Value (Green Channel Maximum), p (8-bit) 


SG2+2YG, © 
4 


/ 


7 


7 


no filter 


SG2+YG, ” 


pixel value max = 255 
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400 mW Blue Laser 


Camera 
Ideal Realistic Quantum 
Photomete 
: 
Reed 
Light 
Meter 


single camera 


| [lux] 


. rY- 





~ dual camera 


—& — . Background (lab lights on) 


— — Laser #4 

---» LaserGlow ($500 laser) 
~--= Laser #1 

— Laser #2 

— Laser #2 with Dust 





10 15 20 25x10" 
time [s 


Simulation of Particle Settling In ee Test Bin 
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These results are generated by using the particle trajectory simulation code (PTQv91), with some minor 
modifications. The plot shows predicted terminal velocity of simulant particles versus diameter. Note that | 
um particles take 10000 s (~ 3 h) to settle | m according to this simulation. That would mean that | um 
particles thrown up to the ceiling (4 m) could take 12 h to totally settle. On the other end of the spectrum, 
The 100 um particles fall one meter in 2 s (8 s from top of the bin). 


v(D) [m/s] 


Dust Particle Terminal Velocity in Regolith Bin 


from PTQv91 Simulation with o = 3100 [kg m3] ea 


10 5 O) 
O vy_m_s_TSE <a 

© vy_m_s__RkK4 CO) O 

—— fit_vy_m_s__RkK4 

— fit_vy_m_s_ TSE _ ee 


: 6 
2 TSE: 
10 z O 
Vp = 7.4131x10°7 7 A 
A = 17.3707 a v(D) = Vo exp (—) 
1+(Do/D) 
Do = 6.077 15 a 
b= 0.483894 Oa 
E 
10 - RK4: 
. : Vy = 1.40929x10-7 
A = 19.6318 
b =0.41719 
107 a 


LASER Side Scatter Model And Measurements 


The following model is a 4-parameter model composed of three sample points along the 
laser beam, as viewed by an oblique angle from the camera. The contour graphs show the 
solution of the 4 parameters - the first 1s the initial intensity value [W/m“%2], and the last 
three are optical extinction coefficients in [1/m]. These values agree well with the laser 
spot measurements. 


fS74000. } 
0.01006 
0.02032 

. 0.0822 | 


emin = 0.0154047 Pp = 





Advantages of this technique as compared to the spot measurement are: 
1) Provides discrete data along the beam path. 
2) Does not require an initial calibration measurement of the no dust case. 


Disadvantages: 
1) More calculations and processing is required. 
2) The results are very dependent on optical alignment (camera relative to laser). 
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Questions 


